Evaluation of heating effects on atoms trapped in an optical trap 






C.W. Gardiner\ J. Ye^'^ , H.C. NagerP and H.J. Kimble^ 

School of Chemical and Physical Sciences, Victoria University, Wellington, New Zealand 
Norman Bridge Laboratory of Physics, California Institute of Technology 12-33, Pasadena CA 91125 
^ Present address: JILA, University of Colorado and 
National Institute of Standards and Technology, Boulder, Colorado 80309-0440 

We solve a stochastic master equation based on tlie tlieory of Savard et al. [T.A. Savard, K.M. O'Hara and 
J.E. Thomas, Phys. Rev. A56, R1095 (1997)] for heating arising from fluctuations in the trapping laser intensity. 
We compare with rec ent experiments of Ye et . al. [J. Ye, D.W. Vernooy and H.J. Kimble, Trapping of single 
atoms in cavity QED, quant-ph/9908 00 7,Phys. Rev. Lett. (1999), in press], and find good agreement with 



the experimental measurements of the distribution of trap occupancy times. The major cause of trap loss arises 
from the broadening of the energy distribution of the trapped atom, rather than the mean heating rate, which is 
a very much smaller effect. 
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In a far-off resonance red-detuned trap, the effective poten- 
tial of the trapped atom can be written 



V{x) = -\a\£ix)\^ 



(1) 



where a is the atomic polarizability, and £{x) is the slowly 
varying field ampHtude [|l|^. Following [J!)], the heating can 
be modeled using a Hamiltonian for a trapped atom of mass 
M of the form 
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which leads to transition probabilities between trap levels of 
the form 
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In these equations, e{t) is a fluctuating quantity, whose spec- 
trum is 



Se{uj) = - dr cos{ujT){e{t)e{t + T)). (4) 
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From these transition probabilities, in follows that the time de- 
pendent probability P{n) that a single atom is in the nth level 
of the trap under the influence of the fluctuation field satisfies 
the stochastic master equation 



P{n) ^^l{n + 2){n + l)P{n + 2) + n{n - l)P{n ~ 2) 
-[n(n-l) + (n + 2)(72+l)]P(n)l. (5) 



with the rate constant 
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FIG. 1. Evolution of the probability distribution P(n): a) Plotted on a short timescale, it can be seen that the heating spreads the initial 
sharp distribution in less than 2ms to cover nearly the full height of the trap; b) Over the full timescale of the experiment losses continue at a 
steady rate. The heating rate used is l/Fe = i/p^xiai ^ 23ms. 



As shown in [y, this constant is equal to the mean heating 
rate, defined as the rate of increase of the level number (pro- 
portional to the energy) of the atom in the trap, i.e.. 



d{n) 
dt 



= r,(n). 



(7) 



It should be noted, however, that this heating rate arises as the 
difference i?„_|_2^„ — Rn-2^n^ in which the quadratic terms 
cancel. If n is significantly different from zero — perhaps 
about 50 in [S] — the positive and negative contributions to the 
heating rate will both be very much larger than the heating rate 
itself. Thus the result of the heating process will be principally 
to spread the distribution over the energy levels, superimposed 
on a much slower increase in the average energy according to 
(0). In fact, the principal time constant for the growth of a, 
the standard deviation of n, is 3Tf/2. 

The principal effect of the heating in the experiment of [||] 
is to expel the atom from the trap, and in general this will oc- 
cur not as a result of the increase of the average energy, but 
rather as a result of the rapid spreading of the width of the dis- 
tribution, so that the upper part spreads into untrapped levels. 

The three-dimensional trap used in [|3| was sinusoidal lon- 
gitudinally, and had a Gaussian form radially. Approximating 
both of these by harmonic fluctuation traps, it was found by 
measuring the fluctuation spectrum that 
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FIG. 2. a) Solid line: Computed probability for the atom 
to remain trapped when the initial mean excitation is the 45th 
level — heating rate as in Fig jit Points: Experimental data from |b|; 
Dashed line: exponential fit to data, b) Solid line: Mean excitation 
of an atom remaining in the trap; Dashed line: Standard deviation of 
the excitation. 

We may safely neglect the much slower radial heating, and 
treat the trap as one dimensional. The trap depth corresponds 
to some 100 levels, so we will model the escape process by 



truncating the master equation to the first 100 levels — once the 
atom leaves this range it ia assumed not to return. The equa- 
tion is easy to solve. As an initial condition, we assume the 
atom is evenly distributed between the levels iVo and A^q + 1 , 
withO < No < 100. The results of a simulation with A^o = 45 
are shown in Fig. 1 . The very rapid spreading of the probabil- 
ity distribution from its initally sharply peaked form is very 
clear In fact very little difference results if a less sharply 
peaked initial distribution is used, even for a width of about 
20 levels. The probability that the atom remains in the trap is 
plotted in Fig. 2a, and this fits the experimental data remark- 
ably well. However, the result is not exponential, though there 
is a strong similarity. Points to note are 

• From Fig. 1 and Fig. 2 it can be seen that a population 
around n = is rapidly produced, and this decays very 
slowly, because the relevant transition probabilities are 
very small. That this is not observed in practice may be 
the result of the existence of other heating mechanisms. 

• The heating rate Fg does correctly give the timescale of 
the heating process, even though the details of the heat- 
ing process are not themselves well summarized by (Q). 

To counter this heating effect one can conceive of introducing 
some kind of laser cooling. One would expect that provided 
the cooling time is sufficiently smaller than the heating time, 
one should be able to ensure that the atom remains trapped. 
We can model cooling by use of a standard master equation 
coupling to a heat bath, such as in [Q], which would give an 
additional contribution to the stochastic master equation (m 
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FIG. 3. Evolution of the probability distribution P{n) 
with both heating and cooling. The heating rate used is 
l/Fc = i/rj"'''' = 23ms, and the coohng rate is Fcooi = 2ms 
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Un+ l)[{n + l)P{n + 1) - nP{n)\ 



+N[nP{n -l)-{n + l)P{n)]\ 



(10) 



In this equation the effective temperature of the heat bath is 
determined by the mean excitation, N, which the bath acting 



by itself would produce in the trap, and Fcooi is the inverse 
cooling time. Adding this cooling term to the heating from 
(g), we see in Fig.^] that the cooling very rapidly counteracts 
the heating. However, in Fig.^ we note that even with quite 
strong cooling, corresponding to l/Fcooi ~ 2ms, the prob- 
ability of remaining in the trap after 60ms is only 90%. By 
solving the equations using only the cooling part (QO), it can 
be verified that most of the loss is in fact a residual effect of 
the heating. 
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FIG. 4. a) Solid line: Computed probability with both heating 
and cooling, for the atom to remain trapped when the initial mean 
excitation is the 45th level — heating and cooling rates as in Fig.H.; 
Points: Experimental data from [H]; Dashed line: exponential fit to 
data, b) Solid line: Mean excitation of an atom remaining in the trap; 
Dashed line: Standard deviation of the excitation when both heating 
and cooling are present. 

However we cannot ensure better trapping simply by in- 
creasing the cooling rate, since the cooling has the effect of 
cooling to a certain residual temperature, and at any non- 
zero temperature there will always be some probability of es- 
caping from the trap, even in the absence of the heating ef- 
fect. Increasing 1/Fcooi at fixed TV (i.e., fixed temperature) 



is equivalent to reducing the timescale of the dynamic pro- 
cesses involved. Once the cooling is fast enough to over- 
whelm the heating, any further increase will simply speed 
up the residual process of trap loss. The only way to get 
more effective confinement is then to reduce the temperature 
to which one cools. With this model of cooling and with 
TV = 10, one finds that the best confinement is obtained with 
1/Fcooi ~ 1ms, although this is only marginally better than 
the case of 1/Fcooi ~ 2ms shown in the figures. 

In conclusion one should bear in mind that the model of a 
truncated harmonic trap is very crude. In the case considered 
here the noise is of the order of 20% of the signal, which also 
means that the validity of the perturnbation theoretic calcula- 
tion used by [11^] to derive the transition probabilities O) will 
also be marginal at best. However the only realistic alterna- 
tives to this very simple picture would involve extensive nu- 
merical work, such as direct simulation of a stochastic differ- 
ential equation, or detailed computations of spectra and matrix 
elements for the appropriate potential. 
Acknowledgements: This work was funded by the Royal 
Society of New Zealand under the Marsden Fund cotract 
PVT902; by the NSF, by DARPA via the Quantum Informa- 
tion and Computing Program administered by ARO, and by 
the ONR. HCN and J Y are supported by a MilUkan Prize Post- 
doctoral Fellowship. 



[1] T.A. Savard, K.M. O'Hara and J.E. Thomas, Phys. Rev. A56, 

R1095 (1997) 
[2] J.D. Miller, R.A. Cline and D.J. Heinzen, Phys. Rev. A47, R4567 

(1993) 
[3] J. Ye, D.W. Ve mooy and H.J. Kimble, T rapping of single atoms 

in cavity QED, |quant-ph/9908007| , Phys. Rev. Lett. (1999), 

in press 
[4] C.W. Gardiner, Quantum Noise (Springer, Berlin 1991); see 

Sect.6.1. 



